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Abstract 



This note will be part of Chapter 6.5 of the ALICE Physics Performance Report 
(PPR), ALICE Physics: Charm and Beauty, where the capabilities of ALICE for 
the detection of open heavy flavour particles will be described. 
We define here the present ALICE baseline for what concerns the heavy flavour 
production cross sections at LHC and the kinematical distributions of the heavy 
quark pairs. We start by qualitatively assessing the Bjorken x regimes accessible 
with charm and beauty measurements at LHC with ALICE (Section^. In Section[2l 
we report the most recent results (and the uncertainties) of the next-to-leading 
order (NLO) pQCD calculations for the cross sections in pp collisions at LHC 
energies. These results are extrapolated to Pb-Pb collisions in Section 13.11 and to 
p-Pb collisions in Section 13.21 taking into account nuclear shadowing and parton 
intrinsic transverse momentum broadening. Heavy quark kinematics as given by the 
NLO pQCD calculation are reported in Section ^ We tuned the PYTHIA event 
generator in order to reproduce such results for what concerns the c and b quarks 
transverse momentum distributions (Section ISJ. Finally, we report the yields and 
transverse momentum distributions for D and B mesons (Sectional). 



1 Accessible x range with heavy quarks at LHC 

In the inelastic collision of a proton (or, more generally, nucleon) with a particle, 
the Bjorken x variable is defined as the fraction of the proton momentum carried by the 
parton that enters the hard scattering process. The distribution of x for a given parton 
type (e.g. gluon, valence quark, sea quark) is called Parton Distribution Function (PDF) 
and it gives the probability to pick up a parton with momentum fraction x from the 
proton. 

The LHC will allow to probe the parton distribution functions of the nucleon and, 
in the case of proton-nucleus and nucleus-nucleus collisions, also their modifications in 
the nucleus, down to unprecedented low values of x. In this paragraph we compare the 
regimes in x corresponding to the production of a cc pair at SPS, RHIC and LHC energies 
and we estimate the x range that can be accessed with ALICE pP for what concerns 
heavy fiavour production. This information is particularly valuable because the charm 
and beauty production cross sections at the LHC are significantly affected by parton 
dynamics in the small-x region, as we will see in the following sections. Therefore, the 
measurement of heavy flavour production may provide information on the nuclear parton 
densities. 

We can consider the simple case of the production of a heavy quark pair, QQ, 
through the leading order^^ gluon-gluon fusion process gg QQ in the collision of two 
ions (Ai,Zi) and (A2,Z2). The x range actually probed depends on the value of the 
centre-of-mass (c.m.s.) energy per nucleon pair ^/s^, on the invariant mass^^ ^qq 
QQ pair produced in the hard scattering and on its rapidity ygq- If the parton intrinsic 
transverse momentum in the nucleon is neglected, the four-momenta of the two incoming 



gluons are (xi, 0, 0, Xi) 



3/2 and {x2, 0, 0, 



X2 are the momentum fractions carried by the gluons, and 



3/2, where Xi and 
is the c.m.s. energy for 
pp collisions (14 TeV at the LHC). The square of the invariant mass of the QQ pair is 
given by: 
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and its longitudinal rapidity in the laboratory is: 
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From these two relations we can derive the dependence of xi and X2 on colliding 
system, Mqq and Vqq: 




A2 



exp - 



(3) 



which simplifles to 



Xi 



exp 



X2 



exp 



(4) 



for a symmetric colliding system (Ai = A2, Zi = Z2). 

Leading order (LO) is 0{al); next-to-leading order (NLO) is 0{a^). 

For two particles with four- momenta {Ei,pi) and {E2,P2), the invariant mass is defined as the modulus 
of the total four-momentum: M = y/ {Ei + E2Y — {pi + p-z)"^- 
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Table 1: Bjorken x values corresponding to charm and beauty production at threshold at 
central rapidity. 



Machine 


SPS 


RHIC 


LHC 


LHC 


System 


Pb-Pb 


Au-Au 


Pb-Pb 


PP 




17 GeV 


200 GeV 


5.5 TeV 


14 TeV 


cc 


X ~ 10-1 


X ~ 10-2 


X ~ 4 ■ 10-^ 


X ~ 2 ■ 10-^ 


bb 






,x ^ 2 ■ 10-3 


X ~ 6 ■ 10-^ 



At central rapidities we have xi ~ X2 and their magnitude is determined by the 
ratio of the pair invariant mass to the c.m.s. energy. For production at threshold [Mcc = 
2 rric — 2.4 GeV, M^-^ = 2 mb — 9 GeV) we obtain what reported in Tabled The x regime 
relevant to charm production at the LHC (~ 10"^) is about 2 orders of magnitude lower 
than at RHIC and 3 orders of magnitude lower than at the SPS. 

Because of its lower mass, charm allows to probe lower x values than beauty. The 
capability to measure charm and beauty particles in the forward rapidity region {y ^ 4) 
would give access to x regimes about 2 orders of magnitude lower, down to x ~ 10"^. 

In Fig. ^we show the regions of the (xi, X2) plane covered for charm and beauty by 
the ALICE acceptance, in Pb-Pb at 5.5 TeV and in pp at 14 TeV. In this plane the points 
with constant invariant mass lie on hyperbolae (xi = M^^/(x2 Snn)), straight lines in the 

log-log scale: we show those corresponding to the production of cc and bb pairs at the 
threshold; the points with constant rapidity lie on straight lines (xi = X2 exp(+2 i/gg)). 
The shadowed regions show the acceptance of the ALICE barrel, covering the pseudora- 
pidity range |?7| < 0.9, and of the muon arm, 2.5 < 77 < 4. 

In the case of asymmetric collisions, e.g. p-Pb and Pb-p'^\ we have a rapidity shift: 
the centre of mass moves with a longitudinal rapidity 

l/c.m. = ^ In f V (5) 



2 VZ2A1 

obtained from equation Q for xi = X2. The rapidity window covered by the experiment 
is consequently shifted by 

Ay 

y\ah. system Vera, system Vc.m.y (6) 

corresponding to +0.47 (—0.47) for p-Pb (Pb-p) collisions. Therefore, running with both 
p-Pb and Pb-p will allow to cover the largest interval in x. The c.m.s. energy in this 
case is 8.8 TeV. Figure |21 shows the acceptances for p-Pb and Pb-p, while in Fig. |21the 
coverages in pp, Pb-Pb, p-Pb and Pb-p are compared for charm (left) and beauty (right). 

These figures are meant to give a first idea of the regimes accessible with ALICE; 
the simple relations for the leading order case were used, the ALICE rapidity acceptance 
cuts were applied to the rapidity of the QQ pair, and not to that of the particles which 
are actually detected. In addition, no minimum pt cuts were accounted for: such cuts 
will increase the minimum accessible value of Mgg, thus increasing also the minimum 
accessible x. These approximations, however, are not too drastic, since there is a very 
strong correlation in rapidity between the initial QQ pair and the heavy flavour particles 

When we write p-Pb, we mean that the proton moves with > 0; when we write Pb-p, we mean 
that the proton moves with pz < 0. 
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M^g = 2.4 GeV: charm □ Central barrel: hi < 0.9 

M™^ = 9 GeV: beauty ■ Muon arm: 2.5 <r| < 4 




"2 



Figure 1: ALICE acceptance in the (xi, X2) plane for heavy flavours in Pb-Pb (left) and 
in pp (right). The figure is explained in detail in the text. 



M™^ = 2.4 GeV: ctiarm □ Central barrel: |ri| < 0.9 

M^'" = 9 GeV: beauty | Muon arm: 2.5 < r| < 4 




Figure 2: ALICE acceptance in the (xi, X2) plane for heavy flavours in p-Pb (left) and in 
Pb-p (right). 
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Figure 3: ALICE acceptance in the (xi, X2) plane for charm (left) and beauty (right) in 
pp, Pb-Pb, p-Pb and Pb-p. 
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it produces and the minimum pt cut will be quite low (lower than the mass of the hadron) 
for most of the channels studied with ALICE. 

The parton distribution functions x f{x, Q^) in the proton, in the CTEQ 4L parame- 
terization are shown in Fig.lH is the virtuality, or QCD scale (in the case of the lead- 
ing order heavy flavour production considered in this paragraph, = Afgg = s xi X2). In 

the figure the value Q'^ = 5 GeV^, corresponding to cc production at threshold, is used. 
The regions in x covered, at central rapidities, at RHIC and LHC are indicated by the 
shaded areas. 

2 Cross sections in nucleon— nucleon collisions 

In this section we present the status of the cross section calculations in nucleon- 
nucleon collisions and their comparison with existing data up to a c.m.s. energy of ~ 
65 GeV. We then report the results for LHC energies. The extrapolation to heavy ion 
collisions is described in the next section. 

The existing data on total charm production cross section in pp and pA collisions'^-* 
up to ISR energies are compared in Fig.Elwith NLO calculations by R. Vogt 0. In Fig.lHl 
a NLO calculation from Ref. [H] is compared to the data in pp collisions from UAl, CDF 
and DO, for which the b quark production cross section integrated for pt > p™™ is given. 
These measurements are taken in the central rapidity region {\y\ < 1.5 for UAl, \y\ < 1 
for CDF and DO). All the calculations have been performed using the following values 
for the heavy quark masses (rric, mb) and for the factorization and renormalization scales 

rrtc = 1.2 GeV ftp = fifi = 2 fiQ (7) 

The pA results were scaled according to the number of binary nuclcon-nuclcon collisions, in order to 
obtain the equivalent cross section in pp. 
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Figure 5: Total charm production cross section from pp and pA measurements compared 
to NLO calculations with MRS D-' (solid), MRST HO (dashed) and MRST LO (dot- 
dashed) parton distributions. 




Figure 6: Comparison with b quark production cross section integrated over pt > p™™ 
from (a) UAl g] and (b) CDF and DO 0. The NLO calculations are with MRS D-' 
(solid) and GRV HO (dashed) parton distributions. 



for charm, and 

mb = 4.75 GeV IJ'F = I^r = /^o 



for beauty; /iq = y (k.q ~^ 'pIq)/'^ + "^Q approximately equal to the transverse mass of 
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the produced heavy quarks. 

For both charm and beauty the theory describes the present data reasonably welL 
The resuhs for LHC energies {^/s = 5.5, 8.8 and 14 TeV) are reported in Table El 
These values are obtained using the NLO pQCD calculation implemented in the program 
by M. Mangano, P. Nason and G. Ridolfi [Zj (HVQMNR) and two sets of parton distri- 
bution functions, MRST HO and CTEQ 5M1 0, which include the small-x HERA 
results. The difference due to the choice of the parton distribution functions is relatively 
small (~ 20-25% at 5.5 TeV, slightly lower at 14 TeV). We chose to use as a baseline the 
average, also reported in the table, of the values obtained with these two sets of PDF. 

Table 2: NLO calculation [2] for the total cc and bb cross sections in pp collisions at 5.5, 
8.8 and 14 TeV, using the MRST HO and CTEQ 5M1 parton distribution functions. 





^pp[mb] 


<[mb] 




5.5 TeV 


8.8 TeV 


14 TeV 


5.5 TeV 


8.8 TeV 


14 TeV 


MRST HO 


5.9 


8.4 


10.3 


0.19 


0.28 


0.46 


CTEQ 5M1 


7.4 


9.6 


12.1 


0.22 


0.31 


0.55 


Average 


6.6 


9.0 


11.2 


0.21 


0.30 


0.51 



The dependence on the PDF set represents only a part of the error on the theoretical 
estimate. An evaluation of the theoretical uncertainties was done by M. Mangano by 
varying the rric {mi,), /iF and /ir parameters and is reported in Table El [ID]- This table 
shows that, at LHC energies, the theoretical uncertainties span a factor ~ 2-3 in the 
total production cross section of both charm and beauty quarks. In the last column of 
the table we report the ratio of the cross section at 5.5 TeV to that at 14 TeV. Despite 
the large spread of the absolute values, the ratio is much less dependent on the choice 
of the parameters; its value is ^ 0.52 for charm and ~ 0.41 for beauty. In Ref. pUj it 
is shown that also the ratios of the pt-differential cross sections are rather independent 
of the parameters choice. This indicates that pQCD can be used to compare the cross 
sections measured in Pb-Pb collisions at ^/s^ = 5.5 TeV to those measured in p-Pb at 
^Snn = 8.8 TeV and in pp at i/i = 14 TeV. 

Yields in proton-proton collisions at y/s = 14 TeV 

Using a proton-proton inelastic cross section a'^^^ = 70 mb at 14 TeV ^Tj and the 
average heavy flavour cross sections in the last row of Table El we calculate the yields for 
the production of QQ pairs as: 

We obtain 0.16 cc pairs and 0.0072 bb pairs per event. 

3 Extrapolation to heavy ion collisions 

In this section we derive the extrapolation of the cross sections and yields to central 
Pb-Pb collisions first, and then to p-Pb collisions. We also point out the different weight 
of the nuclear shadowing effect in the two cases, for charm and beauty production. 
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Table 3: Charm and beauty total cross sections at NLO with different choices of the 
parameters rric ("i^b), A^f and fiR [TU]. In the last column the ratio of the cross sections at 
5.5 TeV and at 14 TeV is reported. 



parameters 


5.5 TeV 


14 TeV 


ratio 5.5 TeV/14 TeV 


TTIq — i.O ore V 

IJ.R = 2/xo Hf = 2/io 


3.7 


7.3 


0.51 


77? — 1 9 HpV 

[mbl = /^o f^F = 2/io 


9.2 


16.7 


0.55 


w — 1 5 OpV 

l-lR = 1^0 I^F = 2/io 


5.4 


10.4 


0.52 


m — ] H OpV 

1 1 — ±.0 vjrc V 
fJ'R = f^O fJ'F = 2/io 


3.4 


6.8 


0.50 


mb = 4.5 vjreV 
fJ'R = 1^0 fJ'F = f^o 


0.20 


0.51 


0.39 


mb = 4:.io (ore V 

fJ'R = f^O fJ'F = fJ'O 


0.17 


0.43 


0.40 


a^l[mh] mb = 5 GeV 

fJ'R. = fJ'O fJ-F = /^O 


0.15 


0.37 


0.41 


mb = 4.75 GeV 
fiR = 0.5/io fiF = 2/io 


0.26 


0.66 


0.39 


mb = 4.75 GeV 

fJ'R = 2/io IJ'F = 0.5/io 


0.088 


0.20 


0.44 



3.1 Nucleus nucleus collisions 

If no nuclear effects are taken into account, a nucleus-nucleus collision can be con- 
sidered, for what concerns the hard cross section, as a superposition of independent 
nucleon-nucleon collisions. Thus, the cross section for hard processes in heavy ion col- 
lisions can be calculated using a simple geometrical extrapolation from pp collisions, i.e. 
assuming that the hard cross section scales from pp to nucleus-nucleus collisions propor- 
tionally to the number of inelastic nucleon-nucleon collisions (binary scaling). 

Nuclear effects — such as nuclear shadowing, broadening of the parton intrinsic 
transverse momentum (kt) in the nucleon, in-medium parton energy loss, as well as possi- 
ble enhancements due to thermal production in the medium — can modify this geometrical 
scaling from pp to nucleus-nucleus collisions. Such effects are, indeed, what we want to 
measure. We chose to include in the simulation only the nuclear shadowing and the broad- 
ening of the intrinsic kt, since they are well established effects. The first effect modifies the 
total hard cross section, while the broadening of the intrinsic kt affects only the kinematic 
distributions of the produced heavy quarks. Nuclear shadowing can be accounted for by 
recalculating the hard cross section in elementary nucleon-nucleon collisions with nuclear- 
modified parton distribution functions and extrapolating to the nucleus-nucleus case. 

The extrapolation, based on the Glauber model fUE], is derived for the collision 
of two generic nuclei with mass numbers A and B, and numerical examples are given for 
the specific case of Pb-Pb reactions at ^/s^ = 5.5 TeV. 

We are interested in the cross section for a sample of events in a given centrality 
range, defined by the trigger settings. The centrality selection can be assumed to corre- 
spond to a cut on the impact parameter b of the collision: < b < be- The sample of 
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X 



2.5 5 7.5 10 12.5 15 17.5 20 

b, [fm] 



Figure 7: Collision geometry in the plane transverse to the beam line (left). Inelastic 
Pb-Pb cross section as a function of the impact parameter cut b < be (right); for clearness, 
here and in Fig. |H1 only the value corresponding to 5% of the total inelastic cross section 
is explicitly indicated. 



events defined by this cut contains a fraction of the total number of inelastic collisions, 
i.e. of the total inelastic cross section, given by 



The definition of the centrality in terms of fraction of the inelastic cross section is more 
appropriate, since the cross section is directly measured, while the impact parameter 
estimation depends on the model used to describe the geometry of the collision. 

In the following, we consider two options for the centrality selection, corresponding 
to 5% and 10% of the total inelastic cross section. The values of be that give these selections 
are 3.5 fm and 5 fm, respectively. 

The inelastic cross section corresponding to a given centrality selection is found 
integrating the interaction probability up to impact parameter b^. 

<b(M = I = 2^ I {l - [1 -^NNTAB(fe)]^^} (11) 

where the value a^^ = 60 mb was used as the nucleon-nucleon inelastic cross section at 
5.5 TeV ^21 5 and the total thickness function Tab 

TAB(fe) = / dhTAis)TB{s-b) (12) 

(vectors defined as in Fig. [71 left) is expressed in terms of the thickness function of the 
nucleus Tj(s) = / dzpi{z,s) for i = A, B, where is the Wood-Saxon nuclear density 
profile ^1] — the thickness function is normalized to unity: / d'^sTi{s) = 1. In Fig. [3 
(right) the inelastic cross section pijl is shown as a function of be- 
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Figure 8: Cross section for a hard process in Pb-Pb collisions relative to the one in 
nucleon-nucleon collisions as a function of the impact parameter cut b < (left). Yield 
of the hard process in Pb-Pb collisions relative to the cross section in nucleon-nucleon 
collisions as a function of the impact parameter cut b < be (right). 



The average number of inelastic collisions for a given impact parameter b is: 

a^^ ■ AB TABib). (13) 

By replacing the inelastic nucleon-nucleon cross section 0"^^ with the elementary cross 
section for a given hard process cipp'''^, we obtain the average number of inelastic collisions 
that yield the considered hard process: 

al^^-ABTMb), (14) 

and the cross section for hard processes for < 6 < b^. 

10 

For minimum-bias collisions (br = +oo), we have: 



<b'(&c) = (tI;'" ■ 271 fjbdb AB Tab(6). (15) 



<B =%TAB. (16) 

The ratio of the hard cross section in nucleus-nucleus collisions, with a centrality 
cut b < be, relative to the cross section in nucleon-nucleon interactions is (see Fig.|Hl left): 

r^^'^(&c) = = 27r / bdb KB TABib). (17) 



The number (yield) of hard processes per triggered event is: 



Arhard/7 \ '^AB ^hard 




where (Fig. |Sl right) 

n(bc) = ^^^^ . 19 

^ J^^bdb {l-[l-aNNTAB(&)]^^} ^ ^ 

For a 5% (10%) centrahty cut in Pb-Pb coUisions, the yield A^ab'^ is obtained by 
muhiplying the elementary cross sections by 26.6 (23.7) mb~^. 

Cross sections and yields in Pb-Pb coUisions at ^s^n = 5.5 TeV 

We used the EKS98 parameterization ^3] of nuclear shadowing, which corresponds 
to a modification of the parton distribution functions of the nucleon in the nucleus, 
ff^i^x, Q^), with respect to the ones of the free nucleon, ff^{x, Q^): 

where the parton distribution functions are given as a function of the momentum fraction 
X carried by the parton inside the nucleon and of the QCD scale and i = q^, Qsca., g 
for valence quarks, sea quarks and gluons. The shadowing factor R^^ for gluons in a 
^''^Pb nucleus is shown in Fig. El The centrality dependence of the shadowing is weak for 
collisions in the considered centrality range (up to 10% of cr"^°') jTH] and it is neglected 
here. 

The reduction of the cross section due to shadowing amounts to about 35% for cc 
pairs, while it amounts only to about 15% for bb pairs, since, as seen in SectionUl beauty 
production corresponds to larger values of x, that are less affected by the shadowing 
suppression. In Section |3] we will show how nuclear shadowing modifies the heavy quark 
kinematical distributions. 

Table lU reports the charm and beauty total cross sections and yields in pp (with 
and without shadowing) and Pb-Pb collisions at ^/s^ = 5.5 TeV, as calculated with the 
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Table 4: Total cross sections and yields for charm and beauty production in pp and Pb-Pb 
collisions at ^/s^ = 5.5 TeV. The effect of shadowing is shown as the ratio Cshad of the 
cross section calculated with and without the modification of the parton distribution 
functions. For the Pb-Pb case we report the two centrality options: 5% and 10% of the 
total inelastic cross section. 





Charm 


Beauty 


agQ [mb] 


w/o shadowing 


6.64 


0.21 


w/ shadowing 


4.32 


0.18 


Cshad 


0.65 


0.84 




5% (t'°°^ 


45.0 


1.79 


10% 


40.1 


1.59 


atQQ 
^^Pb-Pb 


5% 


115 


4.56 


10% 


102 


4.06 



HVQMNR program. The values shown correspond to the average of the results obtained 
with MRST HO and CTEQ 5M1 parton distribution functions. For the Pb-Pb case we 
report the two centrality options: 5% and 10% of the total inelastic cross section. We have 
^Pb%b(5% a''^^V(^Pb^-Pb(10% a'^'') - 1.12. 

3.2 Proton— nucleus collisions 

For the extrapolation to proton-nucleus collisions we use the geometrical Glauber- 
based method already described for the case of nucleus-nucleus collisions. If we consider 
minimum-bias collisions (with no centrality selection), and we use B = 1 and Tb(s) = 6{s) 
for the proton^\ the total cross section for hard processes (fT^ becomes: 

= C ■ 27r / 6 d6 A Ta(6) = A aj^^^ (21) 

^ 

The number of hard processes per minimum-bias pA collision is: 

Arhard \ ^hard / ^inel /'oo\ 

^pA =Acrpp ycTpA- (22) 

Cross sections and yields in p-Pb collisions at ^/s^ = 8.8 TeV 

Using A = 208 and o"p°!pt, = 1.9 barn [T7j, the yield of QQ pairs per minimum-bias 
collision is: _ _ 

= ■ 0. 109 mb"^ (23) 

As for the Pb-Pb case, the effect of nuclear shadowing was accounted for by using 
the EKS98 parameterization ^Hj- Clearly, the effect is lower for p-Pb, since one of the 
colliding nuclei is a proton: the reduction of the cross sections due to nuclear shadowing 
is 20% for charm and 10% for beauty. 

The cross sections and yields for charm and beauty production in pp (with and 
without shadowing) and minimum-bias p-Pb collisions at ^/s^ = 8.8 TeV are reported 
in TableEl The values shown correspond to the average of the results obtained with MRST 
HO and CTEQ 5M1 parton distribution functions. 

A summary of the production yields and of the average magnitude of nuclear shad- 
owing in the three considered colliding systems is presented in Table IHl 

The proton is assumed to be point-like. 
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Table 5: Total cross sections and yields for charm and beauty production in pp and p-Pb 
collisions at ^/s^ = 8.8 TeV. The effect of shadowing is shown as the ratio Cshad of the 
cross section calculated with and without the modification of the parton distribution 
functions. 





Charm 


Beauty 


c^^lmh] w/o shadowing 


9.00 


0.30 


w/ shadowing 


7.16 


0.27 


Cshad 


0.80 


0.90 




1.49 


0.056 


atQQ 


0.78 


0.029 



Table 6: Summary table of the production yields and of the average magnitude of nuclear 
shadowing in pp, p-Pb and Pb-Pb. 





Charm 


Beauty 


system 
centrality 


pp p-Pb Pb-Pb 
min.-bias min.-bias centr. (5%) 
14 TeV 8.8 TeV 5.5 TeV 


pp p-Pb Pb-Pb 
min.-bias min.-bias centr. (5%) 
14 TeV 8.8 TeV 5.5 TeV 


Cshad 


0.16 0.78 115 
1 0.80 0.65 


0.0072 0.029 4.56 
1 0.90 0.84 



4 Heavy quark kinematical distributions 

Figures ITUl and ITD present the transverse momentum and rapidity distributions, ob- 
tained using the NLO pQCD program HVQMNR, for c and b quarks, respectively. The dis- 
tributions for Pb-Pb and p-Pb are normalized to the cross section per nucleon-nucleon col- 
lision. 

We used the CTEQ 4M set of PDF. We verified that the resuhs given by this 
set lie in between the ones obtained with the more recent CTEQ 5 and MRST sets for all 
the relevant kinematical quantities fU]- For the other parameters to values specified in 
Section 121 were used: rric = 1.2 GeV, /xr = /^i? = 2 /^o for charm and mb = 4.75 GeV, = 
A^F = A^o for beauty. Nuclear shadowing is included via the EKS98 parameterization [T^ . 
The parton intrinsic kt is sampled from a Gaussian distribution with mean and a 
(= yj (/ct )) equal to 1, 1.16, 1.30 GeV/c for charm production in pp, p-Pb and Pb-Pb, 
respectively, and equal to 1, 1.60, 2.04 GeV/c for beauty production in pp, p-Pb and 
Pb-Pb, respectively. These values are taken from Ref. The same parameters are used 
also in the calculations shown in the next section. 

In the case of p-Pb events the rapidity distribution in the centre-of-mass frame is 
plotted; the rapidity distribution in the laboratory frame would be shifted by Ay = 0.47. 

We notice that the pt distributions for pp collisions at 5.5, 8.8 and 14 TeV (top-left 
panel) have essentially the same shape. 

The comparison of the pt distributions for pp and Pb-Pb (and for pp and p-Pb) 
at the same centre-of-mass energy shows that nuclear shadowing affects heavy quark 
production only for relatively low transverse momenta {pt < 5-6 GeV/c with EKS98), 
where the QQ pairs are produced by low-x gluons. This is clearly seen in the ratios of 
the distributions, reported in the insets. The value for the upper limit, ~ 5 GeV/c, of 
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Figure 10: Inclusive c quark pt and rapidity distributions obtained from the HVQMNR 
program. The distributions for Pb-Pb and p-Pb are normahzed to the cross sections per 
nucleon-nucleon colhsion and they include the effects of nuclear shadowing and intrinsic 
kt broadening. 
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Figure 11: Inclusive b quark pt and rapidity distributions obtained from the HVQMNR 
program. The distributions for Pb-Pb and p-Pb are normahzed to the cross sections per 
nucleon-nucleon colhsion and they include the effects of nuclear shadowing and intrinsic 
kt broadening. 
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Pair Creation Flavour Excitation Gluon Splitting 




Figure 12: Some of the PYTHIA processes defined as pair creation, flavour excitation and 
gluon splitting. The thick lines correspond to the hard process, the thin ones to the initial- 
or final-state parton shower. 

the pt-region affected by the shadowing in Pb-Pb collisions can be cross-checked with the 
following simple estimate: for the back-to-back production of a cc pair at central rapidity, 
with transverse momenta pi = pi = ^ GeV/c, we have Q — 2j9t = 10 GeV and x ~ 
Q/^Snn = 10/5500 ^ 2 ■ 10"'^; for these values of x and Q, the EKS98 parameterization 
gives R^^ (defined in Eq. ((201)) — 90%. This suppression is already quite small and it is 
partially compensated by the kt broadening. 

A relevant feature of QQ production in p-Pb collisions is a depletion in the forward 
region (where the proton goes) of the rapidity distributions. This effect is due to the 
shadowing which is biased toward forward rapidities, where the smallest x values in the 
Pb nucleus are probed. 

5 Heavy quark production in Monte Carlo event generators 

The program used for the NLO calculations reported in the previous sections is 
not well suited to be included in a simulation, since it is not an event generator and it 
does not provide parton kinematics, but only inclusive distribution. On the other hand, 
widely used event generators, like PYTHIA [TH] and HERWIG ^H]; are exact only at 
leading order, when only the pair production processes, qq — > QQ and gg QQ (see 
Fig. IT^ . are included. Higher-order contributions are included in these generators in the 
parton shower approach (see e.g. Ref. [201 )• This model is not exact at next-to-leading 
order, but it reproduces some aspects of the multiple-parton-emission phenomenon. In the 
following we will concentrate on the PYTHIA event generator; the version we have used is 
PYTHIA 6.150. We have also investigated heavy quark production in HERWIG, observing 
an incorrect behaviour in the final kinematical distributions of both c and b quarks. We, 
therefore, concluded that HERWIG is not suitable for heavy quark simulations at LHC 
energies. More details can be found in Ref. [TUj . 

In PYTHIA, the processes giving rise to contributions above leading order, like (see 
Fig. IT^ Havour excitation, qQ qQ and gQ — > gQ, and gluon splitting, g QQ, 
are calculated using a massless matrix element. As a consequence the cross sections for 
these processes diverge as p^^'"'^ vanishes^^ These divergences are regularized by putting 
a lower cut-off on p^'^'^'^. The value of the minimum p^^'^'^ cut has a large influence on the 
heavy flavour cross section at low pt, a region of prime interest for ALICE physics and 
covered by the ALICE acceptance. Our approach was, therefore, to tune the PYTHIA 
parameters in order to reproduce as well as possible the NLO predictions (HVQMNR). 
We used PYTHIA with the option MSEL=1, that allows to switch on one by one the 

6) phard jg (jgfijie(j 3,8 the transverse momentum of the outgoing quarks in the rest frame of the hard 
interaction. 
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different processes (see the Appendix for more details). The main parameter we tuned 
is the lower p^'^'^'^ hmit. In this procedure we compared the following distributions of the 
bare quarks: 

• inclusive pt and rapidity distributions of the quark (antiquark); 

• mass of the pair: Mqq = ^J{Eq + EqY — {pq + Pq)"^, where Eq = ^ rriQ + pq is the 
quark energy; 

• Pt of the pair, defined as the projection on the plane normal to the beam axis of the 
QQ total momentum; 

• angle A0 between the quark and the antiquark in the plane normal to the beam 
axis. 

In the simulations for Pb-Pb collisions at ^/s^ = 5.5 TeV the parton distribution 
functions used are the CTEQ 4, modified for nuclear shadowing using the EKS98 [T3] 
parameterization. 

Before presenting the results of the tuning of PYTHIA to reproduce the pQCD 
results at NLO, we show that, with the same input parameters, PYTHIA and the pQCD 
calculation in HVQMNR give exactly the same kinematical distributions for the LO pro- 
cess gg — > QQ. The comparison is reported in Fig. for cc production in pp collisions 
at y/s = 5.5 TeV; the normalization is set to the value of the cross section obtained 




Figure 13: Comparison between PYTHIA results (solid histogram) for the LO process 
gg cc, without parton shower, and corresponding HVQMNR prediction (dashed his- 
togram). The centre of mass energy is y/s = 5.5 TeV. 
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Figure 14: Comparison between charm production in Pb-Pb collisions at y/s^ = 5.5 TeV 
in the HVQMNR NLO calculation and in PYTHIA with parameters tuned as described 
in the text. The triangles show the NLO calculation, the solid histogram corresponds to 
the PYTHIA total production. The individual PYTHIA contributions are pair production 
(dashed), flavour excitation (dotted) and gluon splitting (dot-dashed). 




Figure 15: Equivalent of Fig. [T31for beauty production. 
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Figure 16: Comparison between charm production in pp collisions at ^/s = 14 TeV in 
the HVQMNR NLO calculation and in PYTHIA with parameters tuned as described in 
the text. The triangles show the NLO calculation, the solid histogram corresponds to the 
PYTHIA total production. The individual PYTHIA contributions are pair production 
(dashed), flavour excitation (dotted) and gluon splitting (dot-dashed). 




10 20 30 -5 5 




M(bb) [GeV] p,(bb) [GeV/c] A(|)(bb) [rad] 
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Figure 17: Equivalent of Fig. for beauty production. 



for proton-proton without shadowing (first row of Table E} and the PYTHIA results are 
scaled to this number. 

The results of the tuning to pQCD at NLO are shown in Figs. Eland^l where the 
distributions from PYTHIA and the NLO calculation are compared. In this case the overall 
normalization is set to the value of the cross sections obtained for proton-proton with 
shadowing (second row of Table Ej). Despite the fundamental differences between the two 
models, the agreement is relatively good. However, significant discrepancies are present, 
especially in the A0 distribution for cc pairs. 

A similar tuning of the PYTHIA event generator was done also for the production 
of cc and bb pairs in pp collisions at a/s = 14 TeV. The same set of parton distribution 
functions (CTEQ 4) was used, without the modification for nuclear shadowing. Results 
are shown in Figs. El and El The largest difference with the results obtained for the 
Pb-Pb case is a worse description of the rapidity distribution of charm quarks. This 
is due to a feature of the parameterizations of the parton distribution functions: most 
of them, including CTEQ 4, are valid only down to x = 10~^; below this value the 
behaviour depends on the implementation of the specific parameterization but has no 
physical meaning (e.g. for the CTEQ 4 the gluon density g{x) is kept constant at g{10~^)). 
The rapidity range in which the evolution of the parton distribution functions is reliable 
depends on the c.m.s. energy; for charm production at y/s = 5.5 TeV (14 TeV) this range 
is found to be \y\ < 4.3 {\y\ < 3.4), using equation (jU with xi > 10~^ and X2 > 10~^. 
This feature is not present in the latest CTEQ set of PDF, CTEQ 6 [21], which is evolved 
down to X = 10^^. 

The values of the PYTHIA parameters obtained from the tuning are reported in 
the Appendix. 



6 Hadron yields and distributions 

For the hadronization of heavy quarks we use the default Lund string fragmentation 
model [20] included in PYTHIA (JETSET package). The total yield and the rapidity 
density dN/dy in the central region for hadrons with open charm and beauty in Pb-Pb at 
5.5 TeV (5% cr*°* centrality selection), pp at 14 TeV and p-Pb at 8.8 TeV are summarized 
in Tables [71 IHl and ini respectively. The rapidity densities are calculated in —1 < yiab < 1, 
corresponding to —1.47 < yc.m.s. < 0.53 for p-Pb and —0.53 < yc.m.s. < 1-47 for Pb-p. No 
dependence of the relative hadron abundances on the centre-of-mass energy is observed. 

It is interesting to notice the large ratio of the neutral-to-charged D meson yields: 
N{D^)/N(D~^) ~ 3.1. In PYTHIA, charm quarks are assumed to fragment to D (spin 
singlets: J = 0) and D* (spin triplets: J = 1) mesons according to the number of available 
spin states; therefore, iV(D°) : A^(D+) : N{D*^) : A^(D*+) = 1:1:3:3. Then, the 
resonances D* are decayed to D mesons according to the branching ratios. The difference 
between neutral and charged D mesons arises here: due to the slightly larger 4 MeV) 
mass of the D"*", the D*^ decays preferably to D" and the D*° decays exclusively to D°. 
We have (JS]: 

iV(DO) _ iV(DO^i^^^y) + iV(D*+) X BR{D*+ D") + N{D*^) x BR{D*° D") 

iV(D+) ~ iV(D+i ) + iV(D*+) X BR{D*+ D+) + iV(D*o) x fii?(D*o ^ D+) 
H-3 X 0.68 + 3 X 1 



1 + 3 X 0.32 + 3 X 
3.08. 



(24) 
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Table 7: Total yield and average rapidity density for \y\ < 1 for hadrons with charm 
and beauty in Pb-Pb collisions at ^/s^ = 5.5 TeV. The values reported correspond to a 
centrality selection of 5% 0"*°*. 



Particle 


Yield 


(diV/dy)|,,<i 


Particle 


Yield 






68.9 


6.87 




1.86 


0.273 




71.9 


6.83 


B° 


1.79 


0.262 


D+ 


22.4 


2.12 


B+ 


1.82 


0.251 


D- 


22.2 


2.00 


B 


1.83 


0.270 




14.1 


1.30 




0.53 


0.077 




12.7 


1.22 




0.53 


0.082 




9.7 


1.18 




0.36 


0.050 


Ac 


8.2 


0.85 


7^ 


0.31 


0.047 



Table 8: Total yield and average rapidity density for < 1 for hadrons with charm and 
beauty in pp collisions at y/s = 14 TeV. 



Particle 


Yield 


(dAr/dy)|,|<i 


Particle 


Yield 


(dA^/d2/)|,|<i 


DO 


0.0938 


0.0098 


BO 


0.00294 


0.00043 


D° 


0.0970 


0.0098 




0.00283 


0.00041 


D+ 


0.0297 


0.0029 


B+ 


0.00287 


0.00040 




0.0290 


0.0029 


B- 


0.00289 


0.00043 




0.0186 


0.0018 




0.00084 


0.00012 




0.0176 


0.0020 


b: 


0.00084 


0.00013 


A^ 


0.0113 


0.0013 


A^ 


0.00057 


0.00008 


Ac 


0.0110 


0.0013 


a;" 


0.00049 


0.00008 



Table 9: Total yield and average rapidity density for |yiab| < 1 for hadrons with charm 
and beauty in p-Pb collisions at ^/s^ = 8.8 TeV. 



Particle Yield (dA^/dy)|j^j^^|<i 


Particle Yield {dN/dy)\y^^^\^i 


DO + d" 0.926 0.096 
D+ + D- 0.293 0.030 
+ 0.176 0.018 
A+ + A;" 0.118 0.012 


BO + b" 0.0221 0.0030 
B+ + B- 0.0221 0.0030 
B0 + B° 0.0064 0.0009 
A{^ + A;° 0.0041 0.0005 
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p, [GeV/c] [GeV/c] 

Figure 18: Transverse momentum distributions at mid-rapidity for heavy quarks and 
mesons in Pb-Pb at 5.5 TeV. The distributions are normahzed to the same integral 
in order to compare their shapes. 



We chose to use the relative abundances given by PYTHIA, although, experimen- 
tally, the fraction D'^/D"'" was found to be lower than 3. The value measured in e^e~ 
collisions at LEP by the ALEPH Collaboration is ~ 2.4 [22] • This would reduce by about 
6% the expected yield for the D° mesons. 

Figure ITBl presents the transverse momentum distributions at mid-rapidity {\y\ < 1) 
for c quarks and D mesons (left panel) and for b quarks and B mesons (right panel). 



in Pb-Pb at 5.5 TeV. For > and \y\ < 1, we have, on average. 



0.75p^ and 



^ 0.85 The shape of the transverse momentum distributions for D and B mesons 
was fitted to the following expression: 



1 diV 
Pt dpt 




(25) 



The Pt distributions were studied also for pp at 14 TeV and for p-Pb at 8.8 TeV. The 
results of the fits are reported in Table ITUl together with the average pt of D and B mesons 
in the different conditions. The average pt does not depend strongly on the colliding system 
and on the energy in the centre of mass. On the other hand, we remark that (pt) is larger 
by ^ 10% at mid-rapidity than in the forward region (2.5 < y < 4). These two regions 
correspond to the acceptance of the ALICE detector: barrel, |?7| < 0.9, and forward muon 
arm, 2.4 < rj < A. 
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Table 10: Parameters derived from the fit of the pt distributions of D and B mesons to 
the expression and average value of pt for these particles. 







./sZZ TeVl 


»° fGeV/cl 


n 


{■u^ fGeV/cl 




PP 


1 A 




fit; 


i.oO 


U 


p— r D 


fi fi 
o.o 


no 


70 


l.oo 


Ui/labl ^ -LJ 


Pb-Pb 




9 1 9 


9 7S 


1 . 1 




PP 


1 A 


Z. lo 


O.U4 


1.0 ( 


U 


p— r D 


S S 
o.o 


L.LL 


Q 1 1 

o.ll 


1.00 




Pb-Pb 




1 95 


3 1 7 


1 64 




PP 


14 


6.04 


2.88 


4.90 


B 


p-Pb 


8.8 


6.08 


2.90 


4.89 


(IZ/labl < 1) 


Pb-Pb 


5.5 


6.14 


2.93 


4.89 




PP 


14 


6.45 


3.54 


4.24 


B 


p-Pb 


8.8 


6.49 


3.56 


4.24 


(2.5 < yiab < 4) 


Pb-Pb 


5.5 


6.53 


3.59 


4.24 



Appendix: 

PYTHIA parameters used for heavy quark generation at LHC energies 

In Table ^2 we report the complete list of parameters used in the PYTHIA event 
generator fH] in order to reproduce the inclusive pt distribution for the heavy quarks 
given by the HVQMNR program based on NLO calculations by M. Mangano, P. Nason 
and G. Ridolfi [7j. A detailed description of the parameters can be found in Ref. |18j . 

As specified in Section the main parameter we tuned is the lower -p^^"^^ limit: 
the optimal value was found to be 2.1 GeV/c for charm production and 2.75 GeV/c for 
beauty production, both for Pb-Pb collisions at i/Snn = 5.5 TeV and for pp collisions at 
= 14 TeV. Therefore, one can reasonably assume that the same values can be used 
also for p-Pb collisions at ^Js^ = 8.8 TeV. 

The different values for the partonic intrinsic transverse momentum fct in pp, p-Pb 
and Pb-Pb collisions were taken from Ref. 
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Table 11: PYTHIA parameters used for the generation of charm and beauty quarks in 
pp colhsions at 14 TeV, p-Pb colhsions at 8.8 TeV and Pb-Pb colhsions at 5.5 TeV. All 
non-specified parameters are left to PYTHIA 6.150 defaults. 



Description 


Parameter 


Charm 


Beauty 


Process types 


MSEL 


1 


1 


Quark mass [GeV] 


PMAS(4/5,1) 


1.2 


4.75 


Minimum p^"^""^ [GeV/c] 


CKIN(3) 


2.1 


2.75 


CTEQ 4L 


MSTP(51) 


4032 


4032 


Proton PDF 


MSTP(52) 


2 


2 


Switch off 


MSTP(81) 








multiple 


PARP(81) 








interactions 


PARP(82) 








Initial- and final-state 


MSTP(61) 


1 


1 


parton shower on 


MSTP(71) 


1 


1 


2nd Qj.(Jgj^ 


MSTP(2) 


2 


2 


QCD scales 


MSTP(32) 


2 


2 


for hard scattering 


PARP(34) 


1 


1 


and parton shower 


PARP(67) 


1 


1 




PARP(71) 


4 


1 


Intrinsic 








from gaussian distr. with mean 


MSTP(91) 


1 


1 


a [GeV/c] 


PARP(91) 


1.00 (pp) 


1.00 (pp) 






1.16 (p-Pb) 


1.60 (p-Pb) 






1.30 (Pb-Pb) 


2.04 (Pb-Pb) 


upper cut-off (at 5 a) [GeV/c] 


PARP(93) 


5.00 (pp) 


5.00 (pp) 






5.81 (p-Pb) 


8.02 (p-Pb) 






6.52 (Pb-Pb) 


10.17 (Pb-Pb) 
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